Electron pairing in periodic potentials under an external electric field 
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We study the semiclassical dynamics of interacting electrons in a biased crystal lattice. A complex 
dynamical scenario emerges from the interplay between the Coulomb and the external electric fields. 
When the electrons are far apart, the Coulomb potential may be small compared to the external 
potential and the electrons oscillate with effective Bloch frequencies, determined by the local electric 
field. In the opposite case, nearby electrons either separate or form a bound pair, depending on the 
initial energy compared to the band width. The pair due to the Coulomb field is stable even in the 
absence of the external field. 
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Introduction. The dynamics of quantum electrons in 
solids subjected to a uniform electric field is rather non- 
intuitive. According to the semiclassical picture intro- 
duced by Bloch and Zener, noninteracting electrons do 
not accelerate uniformly in real space but oscillate in- 
stead. These coherent oscillations are known as Bloch 
oscillations (BOs). Much after their theoretical pre- 
diction, electronic BOs were observed in semiconductor 
superlattices. ' BOs persist until electrons lose their 
phase coherence through scattering processes. Among 
the various scattering processes that may affect the co- 
herent motion of carriers, electron-electron interactions 
have their own peculiarities. In this regard, Freericks has 
studied the dynamics of conduction electrons and local- 
ized electrons, which do not move but interact with the 
conduction electrons when they are in the same unit cell.'' 
It was shown that BOs are sharply damped and become 
quite irregular in time in this case. Interaction between 
conduction electrons is expected to have less impact since 
all electrons oscillate with the same frequency. Neverthe- 
less, Hubbard-like interactions between particles in the 
same band also induce the irreversible decay of BOs. ' 

Several works have explored the problem of few par- 
ticles in the BO regime. Khomeriki et al. studied the 
dynamics of few interacting bosons in a periodic lattice 
and subjected to a constant force. They found that for 
strong interaction the BO regime re-emerges with frac- 
tional Bloch periods, which are inversely proportional 
to the number of bosons clustered into a bound state. 
The dynamics of two interacting electrons was discussed 
by Claro et al. within the framework of the Hubbard 
Hamiltonian.-' They concluded that electron-electron in- 
teraction induces time-dependent oscillations whose pe- 
riod depends on the strength and range of the coupling 
only. The dynamics of the electron pair without long- 
range interaction also depends on the initial conditions. 
When initially they are far apart, the dynamics is that of 
the single particle BO, as expected. On the contrary, a 
period doubling is found when the two electrons remain 
close, indicating that the pair behaves effectively as a 



composite particle. 

Usually, BOs in the correlated regime are studied with 
contact interaction in the Hubbard Hamiltonian. 
This approximation is well justified for describing two 
electrons interacting by a screened Coulomb potential 
when the screening length is smaller than the lattice spac- 
ing. However, its applicability is questionable when the 
screening length is large. Although long-range interac- 
tions can be implemented in the Hubbard Hamiltonian," 
the resulting equations are complicated, even for two elec- 
trons only. In this paper we consider two interacting elec- 
trons in a lattice subjected to a constant electric field. We 
study their dynamics within the semiclassical framework 
when they interact by the Coulomb potential. We iden- 
tify new oscillation regimes, that were missed in previous 
studies, due to the long-range nature of the Coulomb po- 
tential. In particular, we find that two electrons form a 
bound pair if the energy of the relative motion exceeds 
the half band width, even in the absence of external field. 

Semiclassical approach. The semiclassical dynamics 
of an electron in a periodic lattice is solely parameter- 
ized by its central position r and its central momentum 
hk. Thus, the equations of motion for two independent 
electrons are Kki = —e£ [i = 1,2), where £ is the ap- 
plied electric field. The group velocity of the electrons is 
given by 



Within the tight-binding approximation, the dispersion 
relation of the simple hypercubic lattice is given by 
E{ki) = ~'^JJ2'Ij.=i cos(fci^^a), where a is the lattice con- 
stant, d is the spatial dimension of the lattice, and we 
assume J > hereafter. 

When the Coulomb repulsion between the electrons is 
taken into account, the local electric field is the external 
electric field £ plus the Coulomb field from the other 
electron, which results in the equation of motion 

hk = -e5 -f - , ~ ""^,3 , i^j , (lb) 
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where e is the dielectric constant of the solid. Since 
the interaction term depends only on the relative coor- 
dinate, it is appropriate to make a canonical transforma- 
tion to total and relative coordinates and quasi-momenta. 
Thus, we introduce r ~ ri ~ r2, R — (ri + r2)/2, 
k = (l/2)(fci — fc2) and K = ki+k2. In order to work in 
more convenient dimensionless units, we make the sub- 
stitutions t ^ Jt/h, R ^ R/a, r ^ r/a, K ^ aK, 
k ak to get 



= 2sin(ii:^/2)cos(fc^) 
4 cos (X^/2) sin(fc^) 



, ^ — ^^^^ 
k = -2F 



(2a) 
(2b) 
(2c) 
(2d) 



where the dimensionless magnitudes F = ea£ / J and 
gc = e^/(Jae) have been introduced. The solution of 
Eq. (2c) is trivial and the result can be inserted in Eqs. 
(2a) and (2b) to reduce the number of equations. 

To estimate the range of validity of the equations of 
motion (2), one has to take into account that in real- 
ity each electron is represented by a wave packet with 
a finite width a, which has to be much larger than the 
lattice constant a, but smaller than the separation of the 
wave packets jri — r2|. In terms of the dimensionless 
coordinates, the validity condition is |r| ^ 1. 

Zero external field. Let us consider first the simple 
case when the external field is absent and restrict the 
discussion to one dimension for the time being. Accord- 
ing to Eq. (2c), K — Kq is a constant of motion. Equa- 
tions (2b) and (2d) then form a closed set of equations 
governed by the Hamiltonian of the relative motion 



H, 



9c 



Xq cos k 



(3) 



where Aq = 4cos(A'o/2). The phase-space trajectories of 
the relative motion are given by the contour lines = 
Eq. These trajectories can be classified according to the 
value of the parameter C = _Eo/|Ao|, i.e., the energy of 
relative motion compared to the effective half band width 
Ao = 4| cos(iiro/2)|. As illustrated in Fig. 1(a), there are 
two qualitatively different regimes. If C < 1, trajectories 
are unbounded in x (blue dashed lines in the plot). On the 
contrary, for C > 1 the trajectories are bounded in x (red 
solid line in the plot), thus corresponding to oscillatory 
solutions. The curve defined by C = 1 is the separatrix 
(black dotted line in the plot) between the two regimes. 

Figure 2 shows time evolution of the three trajectories 
highlighted in Fig. 1(a), obtained from the numerical so- 
lution of Eqs. (2d) and (2b) for Kq ^ 0, ko ~ and three 
values of the initial separation, namely gc/xo = 6, 8, 10, 
resulting in C — 0.5, 1.0, 1.5, respectively. In the unbound 
regime C < 1, the relative momentum never reaches the 
edge of the Brillouin zone (BZ), but converges to a value 
smaller than tt. This results in a finite group velocity 
and a ballistic separation of the two particles [see blue 
dashed line in Fig. 2(d)]. This behavior is similar to the 




FIG. 1. (Color online) (a) Phase-space trajectories in ab- 
sence of external field for Ao = 4, given as energy contours 
of Hamiltonian (3) (dark: low energies; light: high energies). 
Trajectories from the unbound regime (C = 0.5, blue dashed 
line), the oscillating regime (C = 1.5, solid red line) and the 
separatrix C — 1 (black dotted line), (b) Phase-space tra- 
jectories for adiabatically varying parameter A. The curve 
for Ao = 4 is the same as the solid red curve in panel (a), 
i.e., Eq = 1.5Ao. As A changes, the energy adjusts itself ac- 
cording to Eq. (6), conserving the phase space- volume un- 
der the curve. For A = the oscillation comes to a halt at 



x/g, = [El - \l) 



-1/2 



0.223. 



dynamics of two electrons in a uniform medium, and no 
signatures of BOs are found, although they move in a 
periodic lattice under an electric Coulomb field. 

The oscillations in the case C > 1 are anharmonic [Fig. 
2(c)] but they are similar to standard BOs, in the sense 
that they are driven by an electric field (in this case, by 
the Coulomb field due to the other electron) and that 
the relative momentum reaches the edges of the BZ [Fig. 
2(a)]. Therefore, the periodic potential and the Coulomb 
repulsion between the two electrons are responsible for 
their pairing. When C approaches unity from above 
(bounded trajectories), both the period and the oscil- 
lation amplitude tend to infinity; the separatrix [black 
dotted line in Figs. 1(a) and 2] corresponds to C = 1. Its 
asymptotic behavior is characterized by fc = it— 5k, where 
5k ^ t~^/'^. Thus, k comes to rest just at the edge of the 
BZ, where the group velocity vanishes. The asymptotic 
behavior of the particle separation is x{t) ~ t"^/^, i.e., 
sub-ballistic [see black dotted line in Fig. 2(d)]. 

Only when the initial conditions satisfy Eq > Xq, i.e., 
the initial energy is too high to be converted completely 
into kinetic energy, the dynamics displays oscillations in 
X (paired electrons) and an unbounded increase of k. 

After having solved the dynamics of the relative mo- 
tion, one solves the equation of motion (2a) for the center 
of mass. The special case Kq — yields X{t) — X{0), so 
the position and momentum of the center of mass remain 
constant, as expected. 
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FIG. 2. (Color online) Time evolution of the relative mo- 
mentum k [panels (a) and (b)] and the particle separation x 
[panels (c) and (d)] in the difTerent regimes of Fig. 1(a) with 
the same color code. Panel (b) shows 5k = n — k{t) for the 
unbound and the separatrix case on a logarithmic scale, the 
gray lines indicating the asymptotic behavior Sk ~ est. and 
5k oc t~^^^, respectively. Panel (d) shows x(t) and and the 
asymptotic behaviors x{t) ~ t and x{t) ~ t^^^ for the un- 
bound and separatrix case, respectively. 



Pairing in the adiabatic regime. We now consider a 
weak external field, such that Fxq is the smallest of 
all energy scales, i.e., Fx{t) <C 1, gc/x{t). A perturba- 
tive approach, however, is not possible because on long 
time scales K grows without bounds and 4cos(i4r/2) = 
4cos(i^i) performs full oscillations. We can, however, 
consider the adiabatic regime, where cos{Ft) varies on a 
much longer time scale than the dynamics of x and k. 
In other words, cos{Ft) can be considered constant dur- 
ing one cycle of x and k. Thus, in the adiabatic limit 
we can safely replace 4cos(Ft) by a constant A. For a 
given value of A, the phase-space trajectories are given 
as equipotential lines of H^. Consequently 



x{k) = 



9c 



E-c + \ cos k 



(4) 



where E^ is the energy of the relative motion. 

We are interested in the paired regime E^ > A (see 
solid red lines in Fig. 1(a) and Fig. 2). As the parame- 
ter A in the Hamiltonian changes slowly, the phase-space 
trajectories are deformed in time. The energy E,. of the 
relative motion is not conserved because of the time de- 
pendence of A. However, the phase-space area enclosed 
by a trajectory is an adiabatic invariant^ '' 

dkx(k) — — = constant . (5) 

^E^ - A2 

Thus, as A deviates from its initial value Ag = 4, the 
energy of the relative motion varies as 



E,{\) = JEi-Xl + X^ 



(6) 



In Fig. 1(b), phase-space trajectories are shown for dif- 
ferent values of A with E^ according to (6), such that the 




FIG. 3. (Color online) Frequency difference Ato = — wi 
obtained from the maxima of the power spectra of Xi{t) and 
X2{t), shown as error bars for gc/xo = 0.1 (red, close to the 
crosses) and Qc/xq ~ 1.0 (blue, close to the diamonds). The 
symbols (crosses and diamonds) show the shift from the time- 
averaged interaction force approximated as Auj ~ 2gc{x~^), 
which deviates from Aa; ~ 2gcXQ^ if the conditions Fxq 2> 
1, Qc/x are violated. 



phase-space area (5) is constant. As A decreases from 
4 to zero, the amplitude of the (anharmonic) oscillation 
vanishes. For negative values, the oscillation is inverted. 
Note that the symmetry (A, k) — (—A, k + n) is clearly 
observed in the plot. Importantly, although the exter- 
nal field F in principle provides a means of getting rid 
of the interaction energy, our results show that this does 
not happen. Instead, the pair remains bound for weak 
external fields. 

Strong field. As a contrast to the adiabatic regime, 
we now consider the regime where the external field F 
is strong, such that Fxq 3> 1, gc/x is the largest energy 
scale. In this regime, the Coulomb force gc/x'^{t) be- 
tween the two electrons is only a small correction to the 
constant force F. Furthermore, the amplitude 2/F of 
the free BO is much smaller than the particle separation 
X, such that, x{t) can be considered as approximately 
constant. Thus, the two particles perform practically in- 
dependent BOs with frequencies uji F — gc/^o and 
UJ2 ~ F + gc/xQ, the frequency difference satisfying 



UJ2 



F 



uji xq _ 2 
9c 



Fxo 



(7) 



as plotted in Fig. 3 (solid line). Figure 3 shows also 
the frequencies obtained from the full integration of Eqs. 
(2) for Ko = ko = 0, gc/xo = 0.1 and gc/xQ = 1.0, 
which show good agreement with Eq. (7), as long as con- 
dition gc/x'^{t) ^ F is well fulfilled. In the opposite case 
of less separated particles, the respective oscillations are 
not independent any more and get distorted. A good 
approximation for the frequency shift is then the aver- 
aged interaction force 2gc{x^'^) (crosses and diamonds in 
Fig. 3). 

Pairing in higher dimensions. Finally, we address the 
phenomenon of pairing due to repulsive interaction in 
higher dimensions for F = 0. Similar to the results shown 
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FIG. 4. (Color online) The pairing transition in d = 3 
dimensions. The initial conditions are fixed by feo = 0, 
^0 ~ 3c/|a;()| — 4d = {0.9, 1.0, 1.1} x 4d, and a random but 
generic orientation of r(0). (a) projection of the momentum 
to the k^-ky plane, the shading indicates the kinetic energy 
in the first BZ. In the case iJo = 1.1 x 4d, fc(t) is shown until 
tj gc = 4. In the other cases, the value Xxmt^oa k is marked 
with a dot. (b)-(d) real-space time evolution of the Cartesian 
components r,i for the different initial energies. 



in Fig. 1(a) and Fig. 2, we numerically integrate the time 
evolution under the three-dimensional generalization of 
Hamiltonian (3) with initial conditions Kq = 0, feo = 0, 
and |ro| such that the energy Eq of relative motion is 
below, equal to, or above the effective half band width 
A := 'iJ2'li=i I cos(it'o,Ai/2)|- Here, we chose Kq = 0, such 
that A = Ad. The orientation of rg is chosen randomly, 
but the ensuing dynamics is generic. The results are 
shown in Fig. 4. In the low-energy regime Eq < A, k 
converges to a point that is different from the BZ corner, 
resulting, again, in a finite group velocity and a ballistic 
separation [see panel (b)]. In the limiting case Eq — A, 
k{t) converges to one of the corners of the BZ, again with 
the asymptotics 6k j — kj —fijii ^ t^^^^ , where the rij are 
odd integers. The real-space asymptotics Xj{t) ^ t^l'^ is 
sub-ballistic [see panel (c)]. When the half band width 
A is not sufficient to absorb the initial energy Eq^ the 
dynamics is quite irregular but the two electrons remain 
paired due to energetic constraints [panels (a) and (d) of 



Fig. 4]. 

If one or several Cartesian components of the initial 
displacement Tq are exactly zero, they remain zero for 
all times. The effective dimension d is reduced and the 
effective half band width A = 4(i as well. However, these 
configurations are unstable: the slightest deviation of 
the initial orientation grows and the system explores the 
whole three-dimensional space at long times. 

Conclusions. A detailed semiclassical analysis of two 
electrons interacting by the Coulomb potential in a bi- 
ased crystal lattice has been presented. The interplay of 
the Coulomb force and the external electric field leads to 
an intricate dynamics that eventually destroys the har- 
monic BOs of independent electrons. Different dynamical 
regimes of the two electrons were observed, depending on 
their initial separation and the magnitude of the external 
field. When the electrons are far apart, the Coulomb po- 
tential may be small as compared to the external poten- 
tial and they oscillates with effective Bloch frequencies, 
corresponding to the local electric field. 

If the external field is weak, then the electrons ei- 
ther separate without oscillations or, when they are suffi- 
ciently close to each other in the beginning, they oscillate 
in the crystal lattice due to the Coulomb field and form 
a bound state. The reason for this pairing to occur is 
the finite band width in the tight-binding lattice. In or- 
der to separate the particles, the initial energy Eq has be 
converted to kinetic energy, which, however, is bounded 
by the half band width A. Thus, the separation of the 
particles is energetically forbidden if Eq > A. Then, in 
dimensions greater than one, k{t) performs a kind of un- 
bounded random walk, resulting in aperiodic dynamics 
of the particle separation x{t) in real space. We have 
shown that in this situation the role of the external po- 
tential is negligible and their dynamics is governed by the 
Coulomb interaction. 
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